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Abstract: We consider the operators with highest anomalous dimension A in the compact 
rank-one sectors su(l|l) and 5u(2) of AA = 4 super Yang-Mills. We study the flow of A 
from weak to strong 't Hooft coupling A by solving (i) the all-loop gauge Bethe Ansatz, 
(ii) the quantum string Bethe Ansatz. The two calculations are carefully compared in 
the strong coupling limit and exhibit different exponents v in the leading order expansion 
A ~ A'^. We find v = 1/2 and = 1/4 for the gauge or string solution. This strong 
coupling discrepancy is not unexpected, and it provides an explicit example where the 
gauge Bethe Ansatz solution cannot be trusted at large A. Instead, the string solution 
perfectly reproduces the Gubser-Klebanov-Polyakov law A = In particular, we 

provide an analytic expression for the integer level n as a function of the U(l) charge in 
both sectors. 
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1. Introduction 

The AdS/CFT correspondence [||, |2|, ^, ^ is a non-trivial map between two integrable 
theories, string theory on AdS^ x and the maximally supersymmetric = 4 super Yang- 
Mills SU(N) gauge theory (SYM). In the planar limit N oo, the string coupling vanishes 
and the correspondence relates a finite superconformal four dimensional theory and free 
string theory on a non-trivial background. Massive string states are predicted to be dual 
to certain composite operators in the gauge theory, with the string spectrum matching the 
gauge anomalous dimensions. In terms of the planar 't Hooft coupling A = ^ym 
duality is of the weak-strong coupling type. Hence, any test of the correspondence must 
exploit some kind of non-perturbative knowledge on at least one of the two sides. 

The large A limit is particularly interesting since the string side can be controlled in 
the supergravity approximation. A quite general prediction is the scaling E ^ 2 ^Jn A^/^ for 
the energy of level n massive string states as A ^ oo |^] . In the gauge theory, a check of this 
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prediction requires the knowledge of the anomalous dimensions of suitable dual composite 
operators in the nonperturbative regime. Such a formidable task is made possible by 
the integrability properties of = 4 SYM. The scaling operators are eigenstates of the 
dilatation operator D that can be identified with the Hamiltonian of integrable (super) 
spin chains in various sectors closed under perturbative renormalization. The integrable 
operator D can be treated by Bethe Ansatz techniques ||6|, ^, ^. In particular, all-loop 
conjectured gauge Bethe Ansatz (GBA) equations are available in the compact su(2), 
5u(l|l) and non-compact sl(2) sectors ||9|, p^ . 

Unfortunately, the GBA equations are only asymptotically exact. For operators with 
classical dimension L, they predict the exact anomalous dimension up to wrapping terms 
appearing at a certain order increasing with L, for instance terms 0{X^) in the 5u(2) 



sector 1 11]. Due to wrapping terms, the GBA equations are not reliable at strong coupling, 
although in some cases they are believed to give the correct leading term. Remarkably, 
in the 5u(2) sector, a local version of the GBA equations has been proposed in the form 
of a Hubbard-like model [^]; it has been conjectured to be free from wrapping problems, 
but the reconciliation of its strong coupling predictions with string theory is far from 
clear H, H, 0]. 



In general, the gauge and string calculations overlap in BMN-like limits |15] where 
L is large. In this case, it is well known that the perturbative comparison in powers of 
A is plagued by the different order of the limits A — > 0, L — > oo on the two sides. For 
instance, the exact three-loop anomalous dimension of two- and three-magnon operators 
in the near-BMN limit [0| exhibits a three-loop discrepancy when compared with the 



leading curvature correction computed in string theory [ [161 , 17, 18, [1^. Similar three-loop 
discrepancies also occur in the expansion around spinning string solutions [Tl| |. 

Along a different route, one can start from the classical string theory (at large A) and 
derive thermodynamical Bethe Ansatz equations at L ^ oo. The discretization of these 
string Bethe Ansatz equations (SBA) have been proposed to compute the leading 1/L 



effects, i.e. one-loop worldsheet quantum corrections [21, y, 10[. The validity of the SBA 



equations at finite L or small A is not guaranteed and indeed they are known to receive 



several kinds of corrections [22, 24, |2^. These corrections have been evaluated for 



various classes of Frolov-Tseytlin spinning string solutions [||, ||, ||, H, |oH|, || . They 
suggest the emergence of an interpolating set of Bethe Ansatz equations working at all A 
and L [^, Q and hopefully solving the three- loop discrepancies. 

In this paper, we take a complementary approach by comparing the GBA and SBA 
equations at strong coupling. Indeed, it is not totally clear to what extent the GBA 
equations are able to predict the correct results in the string regime A ^ cxd as discussed 



for instance in [^, 14, |3^. We attempt to answer this question for the states with highest 
anomalous dimension in the two compact rank-one su(2) and su(l|l) subsectors where we 
are able to solve the GBA and SBA equations at fixed L and generic A. To appreciate the 
special role of the highest states, we briefly summarize some relevant facts about them. 

In the su(2) sector, the highest state is the so-called antiferromagnetic (AF) opera- 



tor |3(:]. It can be defined in the multiplet of operators with fixed classical dimension L. 
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At the perturbative level, it mixes with the other states and its explicit expression is not 



available in closed form at finite L. In the L 
and give the anomalous dimension 



oo limit, the GBA equations can be solved 



A 



su(2) 



lim 

L— >oo L 



1 + 



TT 
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\ 2-K ■ 



(1.1) 



This expression is obtained by taking the L ■ 
to expand at weak coupling and one obtains 



e^ + 1 

oo limit at fixed A. Hence, it is legitimate 
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On the other hand, the strong coupling expansion of Eq. ( |1 . 1| ) is 



lim — 

L— >oo L 



1 r- 3 



(1.3) 



The same leading term is obtained in the Hubbard model formulation of the GBA equa- 
tions 0. As discussed in || this expansion is formal due to the uncontrolled effect of 
wrapping terms. The usual attitude toward this problem is rather optimistic and Eq.(|1.3|) 
with its ~ A^/^ signature is expected to be correct apart from a possible correction in the 
numerical prefactor l/vr^. As a support to the this scenario, it has been proposed to identify 
its dual string state with a suitable string solution in the spirit of similar correspondences 



found in the lowest part of the spectrum [37|. The slow-string solution described in |38| 
exhibits the X^^'^ scaling of Eq. (|l.3| ) although with a different numerical prefactor. This 
quantitative discrepancy has been attributed in [38| to the subtle double limit X,L — > oo. 



However, the Hubbard model formulation |12] suggests that this scenario is not entirely 
satisfactory. There, limit ambiguities and wrapping problems are absent and nevertheless 
the same prediction A/L \/A/7r^ is recovered including the prefactor [l^, |^]. A more 
convincing proof should at least include the analysis of the solution of the SB A equations, 
valid in the strong coupling limit. 

A similar analysis can be attempted for the highest operator in the other compact 
sector su(l|l) [35|. Here the precise form of the operator is known at finite L and simply 
reads Tr{ip^) where ip is the highest weight component of the Weyl spinor in the vector 



1.1) is not known in this sector. In 



multiplet. Unfortunately, a closed formula like Eq. 
the GBA equations for the highest state are studied at weak and strong coupling in the 
L ^ oo limit. The weak coupling expansion turns out to have a finite but rather small 
convergence radius and is not immediately useful to reach the strong coupling regime. On 
the other hand, the strong coupling expansion is ambiguous and depends on assumptions 
about the large A behavior of the Bethe parameters |21]. Two asymptotic solutions for the 



Bethe momenta have been proposed in [35|. The leading term in the anomalous dimension 
is the same in both cases and scales like A^^^. 

To summarize, the present knowledge on the strong coupling behavior of the maximal 
states in the compact sectors of AA = 4 SYM is (we write only the leading term at large A 
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and define L to be the number of Bethe momenta in both sectors) 

? 

(1.4) 

? 

• ; 

3\/3 

CL as L ^ cx), 

where the superscripts gauge/string label the results obtained with the GBA/SBA equa- 
tions. 

In this paper, we pursue the analysis of these states with the aim of filling the gaps 
in the above predictions. We begin with the 5u(l|l) sector which is particularly favorable 
from the technical point of view. At weak coupling we present high-order results for the 
anomalous dimension computed by the GBA equations. We show that a resummation is 
possible by a non-linear acceleration method, the Weniger algorithm. It permits to evaluate 
the anomalous dimension for rather large values of the coupling A. This leads to results 
supporting the A^/^ asymptotic behavior, although with a different coefficient with respect 
to the proposals in 

To investigate further the flow from weak to strong coupling we present a numerical 
solution of the GBA equations. It is worthwile to emphasize that the analysis is quite 
robust and can be extended to very large values of A following in a clean way the evolution 
of Bethe momenta. Our analysis reveals several subtleties involved in the strong coupling 
expansion of the GBA equations. The final result is simple and we are able to compute very 
precisely the leading strong-coupling term of the anomalous dimension. Not suprisingly, 
the agreement with the weak coupling resummation is quite good. At this point, the 
information about the highest state is similar to what is known in the su(2) sector. We 
have accurately computed the weak coupling expansion and the leading asymptotic A^/^ 
term, but we do not know to what extent the GBA are reliable in the strong coupling 
region. We remark that, in the su(l|l) sector, we do not have a slow-string limit solution 
to be identified with the highest state. 

The next obvious step is to analyze the SBA equations in this sector by the same 
methods. Once again the string Bethe Ansatz equations can be integrated numerically. The 
result is very interesting. All Bethe momenta flow to zero at large A with a simple leading 
term ~ A~^/^, and the asymptotic coefficients can be computed numerically. We 
also determine analytically the prefactor in the leading term ~ A^/^ in the anomalous 
dimension. 

The same analysis can be applied to the su(2) sector. The weak coupling resummation 
does not apply here because the precise form of the highest AF state depends on L. How- 
ever, the numerical and analytical study of the strong-coupling behavior of GBA and SBA 
equations can be performed without difficulty. We find a pattern similar to the su(l|l) 
one. Our results are summarized in the following table which has to be compared with 
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A gauge AString 

lim ^ = ^ AV^ (exact) lim 

^gauge AString 

=CiAV2, (conjecture m) lim 



Eqs. ( |l.4| ). In the right hand sides, we only report the leading term at large A and finite L 



A gauge , . 

1 A string 

2L ^ A , y_L_ _ _ 



27rL sin ^ 2L 2 

^string 



(1.5) 



CL ^ 0.1405(1), asL^oo. 



The ratio A^*"^!^^ / {2L) is independent on L. Also, in the L — > oo limit, we recover the 



exact prefactor l/vr^ for A™/(2L). 



su(2) 

'su(2) 

The physical contents of Eqs. (|1.5| ) will be discussed in Sec. ^, after heving illustrated 
the technical details of the derivation. These will be organized as follows. In Sec. ^, we 
present the gauge Bethe Ansatz equations for the su(l|l) sector. In Sec. ^]l| and |2]^ we 
analyze them for the highest state obtaining in particular the strong coupling expansion of 
the anomalous dimension. Sec. ^presents the string Bethe Ansatz equations, again for the 



su(l|l) sector and Sec. 3.1 and 3.2 repeat the previous analysis. In this case, the strong 
coupling expansion of the anomalous dimension is determined exactly. In Sec. ^ and |5| we 
present a similar analysis for the other compact su(2) sector. 

2. The gauge Bethe Ansatz equations for the highest state in the su(l|l) 
sector 



The dilatation operator in the su(l|l) sector can be associated with a super spin chain |39 
19]. The all-loop gauge Bethe Ansatz equations have been proposed in |lO| and read 

1 



_2x+{pk)x (pj) 

jy^k 1 



,^Lp^ = r] — ^^^^ k = i,...,L, (2.1) 



2x ipk)x+{pj) 
where L G 2N + 1 and 



P 

The coupling g is related to the 't Hooft coupling by A = 8 vr^ (7^ . The anomalous dimension 
of the state associated with the solution {pkid'^)} is 

L 



(2.3) 



The Bethe momenta of the highest state at (7 = are 



27r L - 1 L - 1 , ^ 

Pk = jjnk, nk = 2~''"'^~' 
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The GBA equations in logarithmic forms are 



Pk = ^ nk - - l^log (2.5) 

j^k 1 



where are given in (|2.4| ). They are suitable for weak coupling expansions since the 
second term in the r.h.s. is 0{g^). 



2.1 Weak coupling expansion and Weniger resummation 



The weak-coupling expansion of the anomalous dimension is easily obtained. We simply 
start with the zero-th order value of Bethe momenta for a certain fixed L, Eq. (p.4|). 
Then, we replace them in the GBA equation and expand the r.h.s. at first order in . 
Repeating this procedure and expanding the expression for A order by order, we obtain 
the perturbative expansion of A. 

In the ratio A/L, the terms up to 0{g'^^~'^) do not change if L is increased. For 
instance, 



A\ 3 9 . 4 . fi 235 o 2209 ^ / s 

-J =- + 252.4/ + 14/ __/ + ^5io + ... ^ (2.6) 

A\ 3 2 . 4 29 259 o 2611 ^ / ^^ 

r),., = 2+'»-'"'+T''--« +—"+■■■• 

A\ 3 2 . 4 29 259 o 1307 ^ / ^ 

^j^^^= 2+2,^-4/ + -/--/ + ^,- + ..., (2.8) 



This remark allows one to compute the L ^ 00 limit of the expansion at a fixed order in 
g by simply taking a sufficiently large L. The procedure can be performed at the semi- 
numerical level. In other words, we work with finite high precision using numerical values 
of the momenta in a symbolic algebra calculation. If the precision is suitably high, the 
identification of the coefficients of the energy expansion can be unambiguously identified 
with rational numbers. We have performed the calculation up to the term g^^. This 
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requires L > 35. The result is 



lim ^ = 3+252_454 + ^g6_ 259 8 +1^ 10 _ 1790^12 + 103965^4 _ 5(^97 16 
6324557 „18 _ 40702709 „20 , 8561442701 „22 _ 114529021311 „24 

16 y 16 y 512 y 1024 y 

777307887947 26 _ 2670717561365 28 , 37098574647961 30 _ 4161069724993527 32 
1024 y 512 y 1024 ^ 16384 ^ 



+ 
+ 



' 1024 y 512 y ' 1024 y 16384 » 

I 29408079892945107 34 _ 104670245742870895 36 , 5999052730939686071 38 
16384 y 8192 ^ 65536 ^ 

86452214868942845981 40 , 626162974135003430373 42 _ 4556537471418865642837 ' 
131072 y 131072 ^ 131072 ^ 

I 66598702591887298874029 „46 _ 244303906058015917431755 „48 

262144 y 131072 ^ 

I 7195137546781961772111605 „50 _ 106303776929607820974312023 _52 

524288 y 1048576 ^ 

I 49230031886653815687152661 54 _ 45726319914455572899079305 56 

65536 y 8192 y 

I 348843198908576206971428650203 58 _ 2605279742772089252587976183821 60 

8388608 y 8388608 y 

I 39003030225010830621366145740085 62 _ 2340608578131628813286501122058923 J 
16777216 y 134217728 y 



16777216 y 
I 35185861176795745832756768610959237 66 
268435456 y 



268435456 • 
264968465576189708105542064159612145 
268435456 



+ (2.9) 



We have identified the rational numbers by working with 200 digits arithmetics ^. This 
power series is convergent for 5^ ^ |, a rather small convergence radius. It is definitely 
useless to evaluate strong coupling behavior at least in this form. We need an analyti- 
cal continuation beyond the convergence radius. Since the series is alternating, we have 
attempted such continuation by means of the non- linear Weniger algorithm |41| that we 
describe in Appendix |^. 

In our case, the Weniger algorithm is found to work very well. As an example, we 
consider = 1 which is far beyond the convergence radius. The Weniger approximants 
are shown in the first curve of Fig. (|^) where a clear and definite convergence is achieved. 
Going to higher values of the coupling, we find that the resummation algorithm gives stable 
results up to g"^ ~ 20, which is a fairly high value. The other two curves of Fig. (||) shows 
the behavior of Weniger approximants for g^ = 10, 30. 

The plot of A/L — 3/2 at L ^ oo in the stability region is shown in Fig. (^) where we 
also draw the results from the numerical analysis of the GBA equations that are discussed 
in the next Section. From the Weniger algorithm, we recover clearly the asymptotic \/A 
behavior. A numerical fit gives the estimate 

lim — ~ 0.1404 \/A. (2.11) 



^For instance, with such a precision, the coefficient of g^^ appears in the calculation as a floating point 
number r such that 

268435456 r = 264968465576189708105542064159612145. 000^ ■ ■ , (2.10) 

163 null digits 

allowing for a safe identification of the numerator of the rational coefficient. Standard packages like 
Mathematica |t40| allows easily this kind of high precision numerical calculations. 
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The numerical prefactor is different tlian the one predicted in |35]. The question is whether 
the resummation is failing or the strong coupling expansion is revealing some surprise. In 
the next Section, we shall answer this question in favor of the second hypothesis. 

2.2 Numerical solution and strong coupling behavior 
2.2.1 Preliminary remarks 

The form of the GBA equations at strong coupling depends crucially on certain a priori 
assumptions about the asymptotic form of Bethe momenta. For the highest state, it is 
necessary to consider separately three different cases, motivated by the following numerical 
analysis. If we denote by p{g) a particular running Bethe momentum, we consider the 
three special large-5 behaviors 



I 

II 
III 



pig) 
p{g) 
p{g) 



>p>0, 

ag~'l\ 



(2.12) 



ag 



-1 



The ratio g/x {p) has the following limit 



I: V2e^' 



2 e(sm§j, 



x^{p) 



II: ^/2e{a), 
111: 2" 



(2.13) 



1 + Vl + 2a2' 



where e(x) = x/\x\ is the sign function. For better uniformity, it is convenient to write the 
case II as _ 



V2e{a) = \/2e^*2 e(sin|) 



(2.14) 



where p = in this case, and e(0) = ±1 according to the sign of a. With these conventions, 
cases I and II are expressed by the same formula. 

2.2.2 The Arutyunov-Tseytlin Ansatz 

An Ansatz for the strong coupling behavior of the Bethe momenta of the highest state 
is described in |Q. As we shall discuss, it is closely related to the actual solution. The 
Arutyunov-Tseytlin Ansatz assumes that one p remains zero and the other tend to non- 
zero limits symmetrically distributed around zero. This Ansatz is quite reasonable since 
the symmetric pattern is valid at 17 = and remains true at all orders in the weak coupling 
expansion. 

Under this assumption, we can look at the positive p only. Using the expressions (|2.13| ), 
the GBA equation for any of them reduces at strong coupling to 



iLpk 



L-l 
2 

n 



1 _ f,-i(Pk-Pj) I ^ g-5(Pfe+Pj) 



62 



kiPk-Pj) 



1+62 



HPk+Pj) 



2 

n( 



(2.15) 
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This gives 

L-3 L- 1 , , 

Lpk= ^ TT —pk + 27rmk, mfe G Z, (2.16) 

or 

Pi- = h 7r . 2.17 

^ 3L - 1 3L - 1 ^ ^ 

2.2.3 Explicit results at various L 



As we discussed, the result ( 2.17 ) implies an asymptotic anomalous dimension which does 
not agree with the resummation results. To understand what is happening, we have solved 
numerically the gauge Bethe Ansatz equations according to the following recipe 

1. we start with the solution at g = 0, 

2. we progressively increase g and solve step by step the Bethe equations by Newton's 
algorithm |42], 

3. at each step, we use the solution at the previous 5 as a starting guess for Newton's 
algorithm. 

The procedure turns out to be very stable and can be extended up to very large A values. 
In particular it is possible to increase g in logarithmic scale. The stability of the algorithm 
is checked by varying the numerical precision used in the intermediate computations. We 
never encountered any singularity. By means of this numerical method, we have investi- 
gated the GBA equations at several L in order to discover why and when the Ansatz ( 2.17] ) 
fails. 



Data for L = 3 and L = 5 are reported in Fig. (^); they confirm the Ansatz ( 2.17 ). We 
show the positive Bethe momentum and the prediction with k = 1 for L = 3, and k = 0,2 
for L = 5. Notice that convergence is achieved at quite large A. 

At L = 7, the Bethe momenta are shown in Fig. Here, something new happens. 
One of the Bethe momenta tends to zero. Nevertheless, the Ansatz ( ^.17 ) is still working, 



with k = 0,2. The reason is that the vanishing momentum tends to zero like A~^/^ (case 
II) and the limiting form of the GBA equations is the same as it would be in case I. The 
asymptotic form of the vanishing momentum is illustrated in Fig. (|5|). 

At L = 9, the Bethe momenta are shown in Fig. (|6|). Here again one of the momenta 
tends to zero. Now, the Ansatz (2.17) fails to predict the correct asymptotic values of the 



non- vanishing momenta. Indeed, the vanishing momentum tends to zero like A~^/'^ (case 
III) as shown in Fig. (0) and the limiting form of the GBA equations is changed. The 
dashed lines predicting the actual asymptotic non-zero p are obtained as follows. Using 
again Eqs. (|2.13| ), the GBA equation for any positive p reads in the strong coupling limit 
(L = 9) 

where a appears in the asymptotic form of the vanishing momentum which is a /VX. For 
each p we can determine the three positive p nearest to the numerical asymptotic values. 
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Then, we fix the parameter p by using the strong couphng limit of the GBA equation for 
the vanishing momentum 



i=i,2,3l-/oe *2 l + pe*2 
The numerical solution is easily found. With 25 digits, it reads 

p = 0.7261948032180057677773276, (2.20) 
pi = 0.6047720145641787805731663, 
P2 = 1.648738996485669279031225, 
P2. = 2.646257150202204974776960. 

The agreement, shown in Fig. @ is excellent. In the following of this paper, we shall often 
have to solve equations like Eqs. ( p.l8 - 2!T9|) . Whenever we claim that a numerical solution 



is easily found, we mean that standard packages, like Mathematica |40|], can determine the 
solution with high precision in a straightforward way. For simplicity, we shall give just a 
relatively small number of digits for such results, but in all cases, we have checked their 
stability by increasing the precision and checking that the result is unchanged. 

At L = 11, the Bethe momenta are shown in Fig. (P). Here again one of the momenta 
tends to zero like A~^/^ (case III) and the limiting form of the GBA equations is changed. 
As before, we can compute the dashed lines showing the asymptotic non-zero p. The GBA 
equations for the positive p read (L = 11) 



1 _ /)2 p-ip 

eiLp = -e-^^P- — , (2.21) 
1- p^e'P ^ ^ 



Again, for each p we can determine the four positive p nearest to the numerical asymptotic 
values. Then, p is fixed by the GBA equation for the vanishing momentum which reads 

■El _ -El 

nl — /9e*2l + pe*2 
, ' ' (2.22) 

j=i,2,3,4 1 - pe * 2 1 + pe* 2 

The solution is now 

p = 0.616048, (2.23) 
pi = 0.227432, 
P2 = 1.09891, 
P3 = 1.92557, 
P4 = 2.73741. 

The agreement is shown in Fig. 

At L = 13, the Bethe momenta are shown in Fig. (^. Here two momenta tend to 
zero, one like A~^/^ and the other like A~^/^. We repeat the exercise of computing the 
asymptotic p. The GBA equation for the positive p reads {L = 13) 

1 _ -ip 

^iLp^^-5ip ^L^, (2.24) 
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For each p we can determine the four positive p nearest to the numerical asymptotic values. 
Then, we fix the parameter p as before by the GBA equation for the vanishing momentum 
which reads 



1 



n ^ 

=1,2,3,4 J- 



1 — pe*2 i-|-pe~ 



pe 



l + pe' 



(2.25) 



The solution is 



p 


= 0.53244, 


Pi 


= 0.72414, 


P2 


= 1.42787, 


P3 


= 2.11751, 


Pi 


= 2.80082. 



(2.26) 



The agreement is shown in Fig. (|^. 

If L is further increased, the pattern of vanishing and non vanishing momenta turns 
out to be quite regular. In the following table we show for each L the number of 
positive momenta vanishing like X~^^^ and the number of those vanishing like X~^^'^. 
The general formulas expressing the Z3 regularity of the Table are 



L 

Ni 
2 

±. 



7 9 

1 

1 



11 
1 





13 15 17 
1 2 2 

1 



19 21 23 
2 3 3 
10 



25 27 29 
3 4 4 

10 



Table 1: Periodicity of the number of vanishing Bethe momenta. 



2 



Ni 

4 



L-1 



L- 1 
0, — - — mod 3 = 0, 



1, otherwise 



L-1 
1, — - — mod 3 = 0, 



0, otherwise 



(2.27) 



(2.28) 



For instance, if L = 43 we expect 6 momenta vanishing like A~^/^ and one like X~^^^. The 
full set of momenta is shown in Fig. (p!o[). It can be checked that the vanishing momenta 
have precisely these asymptotic behaviors. 



2.2.4 Asymptotic form of A 

The asymptotic form of the anomalous dimension at fixed L and large A is 

A A- 1 -J^ I pk 

— clVX with CL = -f— > , sm — , 

/, IjTi ^ — ^ I J 



k=l 



(2.29) 
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where the pk are the asymptotic non-zero values of the Bethe momenta. Due to the above 
periodicity, we can estimate Cqo by considering separately our data for cl for the three 
values of ((L — l)/2)mod3. The result from a fit of the data at L > 5 by using a cubic 
polynomial in 1/L are shown in Fig. (11). The three subsequences have clearly the same 



limit. We find 

Coo = 0.1405(1), (2.30) 

where the error is a conservative estimate of the finite L fit. 

A remark about the order of the two limits L, A ^ oo is in order. We applied the 
resummation algorithm to estimate the leading term at large A of lim^^^oo A/L. Here, 
solving the GBA equations, we are fixing L and taking the large A leading term ~ ci,\f\. 
Then, we evaluate Coo = limi^oo c/,. Therefore, the double limit L, A — > oo is taken in two 
different orders. Nevertheless, the agreement of the leading term in the two calculations is 
not surprising. This is precisely what happens for the AF state in the su(2) sector. There, 
one can start from Eq. (^]5) and take after the A ^ oo limit. Alternatively, one can take 
the large A limit at fixed L, e.g. in the Hubbard model formulation. The result for the 



first three terms in the expansion is the same as discussed in |14| 



3. The string Bethe Ansatz equations for the highest state in the su(l|l) 
sector 

The analysis of the GBA equations is certainly interesting, but the ultimate goal is the 
comparison with string theory. As discussed in the Introduction, it is not clear to what 
extent the GBA predicts correct results at strong coupling. This question can be investi- 
gated by studying the string Bethe Ansatz equations |9|, |l^] expected to predict the correct 
strong coupling behavior of string states, at least at large L. 

In order to write the SB A equations in a compact way, we define = x^(pfc) and 
"Wfc = u{pk) where 



u{p) = - cot-Jl + 8g'^ ^^'^^ 2' 



(3.1) 



The string Bethe Ansatz equations are then 



n 



i^k 1 



O ^ 1 



2Xf^ X- 



(3.2) 



where the scattering phase •& is 



9' \ 



-2 



(1 



9' \ 



2i{uk-Uj) 



2X+ X. 



2^fc 



2^fc 4 



9 



2xtx-} J 



(3.3) 
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In logarithmic form we have 



27r i 



9^ 



log + 2 i (-Ufc - Uj ) log ■ 



1 r 1 1 



(3.4) 

These equations are considerably more involved than the GBA ones. Nevertheless, we have 
been able to repeat step by step the previous analysis as we now discuss. 

3.1 Weak coupling expansion and Weniger resummation 

We repeat the semi-numerical algorithm that we followed for the weak coupling expansion 
of the GBA equations. We take again L = 35 and obtain the result 

A string 

i;™ _ 3 1 9 „2 _ ^ „4 I 25 „6 _ 601 „8 I 2849 ^10 _ 25141 „12 , 429809 „14 _ 9022721 ^16 

/ml — ^9^ 2 9 12 9 ^ 12 9 20 y ^ 60 y 210 9 

I 149821573 18 _ 8640293477 „20 , 1812303079883 „22 _ 88730558092937 „24 

560 9 5040 9 161280 9 1182720 9 

I 1804497110708207 „26 _ 26846650998855167 _28 , 279571052498891591 „30 
3548160 9 7687680 9 "t 11531520 9 

125402745492098095339 32 , 9738744677918359729 34 _ 6703245537745284313789 36 
738017280 9 "t 8110080 9 784143360 9 

I 1024647660942740023729097 38 _ 2525524043347946614344579101 40 

16728391680 ^ 5721109954560 ^ 

I 18296543439438265154466562553 42 _ 1331791971895366969823043659509 44 

5721109954560 9 57211099545600 9 

I 9735683703553399316563746438049 46 _ 345964047352595964568301174322841 48 

57211099545600 9 277022166220800 9 

I 48414229366853228412090322467734657 ^50 _ 2689943057205158649159985079500991 _52 

5263421158195200 9 39574595174400 9 

I 5303581496963536951888629403691895353 54 
10526842316390400 9 
19710992154318324101213041165671461761 56 
5263421158195200 9 
I 870504639734489849012829744121725553 58 
31190643900416 9 
10693434560671288576899374290035763214898727 „60 
51286775765454028800 9 
I 76039980287752685716954083401449551327779 62 
48705390090649600 9 
10643182102443547953784209313106048717253607517 64 
908508599273757081600 » 
I 2240647681205846470476844290367422504232816365029 66 
25438240779665198284800 9 
5626138712126728829417020286414076511813854415733 68 i (n(n'J^\ (r>rC\ 
8479413593221732761600 9 ^ '^\9 )■ l-J-^j 

The Weniger resummation algorithm is convergent for < 10. We show the resummed 



expression for A in the left panel of Fig. (12) for the gauge and string cases. The right 
panel shows the derivative 

d ( A^t"°s\ , , 

log hm — - — , (3.6) 



dlogA \L^oo L 

which estimates at large A the exponent of the leading term. In the gauge case, it ap- 
proaches the value 1/2 at large A, as we discussed. In the string case, the asymptotic 
value appears to be definitely smaller and the figure is qualitatively compatible with an 
asymptotic behavior ~ A^/^. 
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To support this conclusion, we have fitted the whole data with the functional form 

A string)' \\ 

lini -A:^ = CO A'' + ci + C2 A-^ + C3 A'^^ (3.7) 

The standard is a quantitative measure of the deviation from the supposed dependence 
on A. The best values for the exponent v are 

gauge BA : v^t = 0.496, string BA : v^t = 0.24, (3.8) 

which are very close to 1/2 and 1/4. If we now fix the exponent = 1/2 or 1/4, we find 
the following for the two curves: 



fixed exponent u = 1/2 v = l/A 



x'Cgauge BA) 2.3 • 10"^ 0.02 ^3 
X^(stringBA) 0.03 1.4-10"^ 



As a conclusion, the resummed anomalous dimension favors the choice = 1/4 in the 
string case. The leading term with its numerical prefactor is 

A string 

lim — — = 0.70(1) A^/^ (3.10) 

L— >00 Li 

The prefactor is difficult to estimate and a better determination would require a stable 
resummation at larger A. 

While this result is quite pleasing, it must be criticized because of the moderate re- 
summation range. Based on the weak coupling arguments [presented so far, it can not be 
ecluded that the string curve in the right panel of Fig. ( p!2|) could rise at larger A and flow 
back to the gauge value 1/2. To pursue the analysis, as in the gauge case, we turn to a 
numerical iterative solution of the SBA equations. Indeed, it should be clear that several 
solutions are possible at strong coupling and the problem is again that of choosing the right 
one. 

3.2 Exact solution at strong coupling 

We now determine the numerical solution of the SBA equations. We follow the same 
procedure we described for the gauge BA equations. The result is fully consistent with the 
weak coupling resummation: All Bethe momenta pk vanish at large A with an asymptotic 
behavior pk ~ Ofc A~^/^ for all pk\ We illustrate this noticeable result by showing in Fig. (13) 
the evolution of Bethe momenta scaled by A^/^ in the four cases L = 3, 5, 15, 29. 

The coefficients are symmetrically distributed around zero. Qualitatively, L — 2 
coefficients are almost evenly spaced around zero. Two special momenta have instead 
coefficients well separated from the central band. Looking in more details at the explicit 
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solution, one finds that the in the central band have a non-trivial asymptotic density. 
The analysis of the large L form of this density if deferred to future work. 

It is not difficult to find an exact equation for the asymptotic coefficients a^. As we 
said, one of them is zero, (L — l)/2 are positive, and {L — l)/2 are opposite to the positive 
ones. Expanding the SBA equations at large A is a bit tricky but straightforward. We 
obtain the following equation determining the positive > 



exp I 7^- Ofc ai 



2tt 



n 



«i) + 2vr(afcVa.i) 



aj>0 2 

Ik 



^{ak - aj) + 27r(Q^ + a 



aj) + 27r(a^ ^ + a,-^) 



4j(/i(ofe)-h(aj)) 



(3.11) 



h{a) = — - — . 

Ibvr 



where in the l.h.s. the sum includes the case j = k. 

Although Eq. ( 3.11| ) is rather complicated, it can be solved numerically without dif- 
ficulties, at least starting from the numerical p obtained at a reasonably large A. As an 
example, at L = 5, 7, 9 we find the following (numerical) solutions 



L 



L = 7 



L = 9 



02 



2.9213116645, 
7.9614845209, 



ai 
as 



2.1234902933, 
4.0857786417, 
9.1813290270, 



ai 

04 



1.6905819725, 
3.1495210704, 
4.8259074880, 
10.203166001 



(3.12) 



These values are in perfect agreement with the numerical solution of the string Bethe 
Ansatz equations as illustrated in Fig. (14). 

Actually, if we are interested in the asymptotic form of A''*''™^, we do not need the full 
information encoded in the a^, but just their sum. Indeed, 



. string 



L 



-y 



(3.13) 



where 

1 ^ 

^^ = ^Y.\^^\ = L^Y.^^- (3-14) 

k=l afe>0 

We now take the product of Eqs. ( |3.11| ). The right hand sides cancel perfectly. Evaluating 
the product of the left hand sides we obtain 



exp 



(2^)2 iVi, 



(3.15) 



where Nl G N. These integers can be determined by solving Eqs. ( 3.11| ) at a certain L. The 
starting point for Newton's algorithm is taken from the solution of the SBA equations at a 
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reasonably large A. Given the solution for the coefficients, we can compute N^. Indeed, 
the solution of Eqs. ( 3.11| ) can be accomplished easily with an arbitrarily high number of 



digits and the identification of the integer Nl is totally straightforward and unambiguous. 
The first values of A'^^;, are 



L : 5 7 9 11 13 15 17 19 21 
TVl : 3 6 10 15 21 28 36 45 55 



(3.16) 



The following simple formula holds 



leading to the prediction 



Nl = 1{L^-1), (3.17) 



1 1 1 / 1 \ "'^''^ 

Now, we do not observe any particular Z3 structure. The values of cl are perfectly smooth 
as L increases. As a further check of this analytical expression, we show in Fig. (^) the 
fit of Coo with a simple quadratic polynomial in 1/L. There is perfect agreement with the 
prediction 

Coo = (3.19) 

The numerical solution of the SBA equations and the resummation in the stable region 
9^ ^ 10 in perfect agreement. 

Of course, the appearance of the integer in the asymptotic form of A at large A is 
not surprising. Indeed, this is a general feature of the SBA equations as discussed in pl| . 
If the Bethe momenta vanish like A~^/^ at large A, then the asymptotic form of A is 

A~2V^Al/^ (3.20) 

where n is a sum of mode numbers. This is the celebrated string prediction of where n 
is the level of a massive string state. The calculation that we have described identifies the 
precise value of n = Nl for the state dual to the highest operator in the 5u(l|l) sector. 

4. The gauge Bethe Ansatz equations for the AF state in the su(2) sector 

It is straightforward to extend the analysis to the highest state in the 5u(2) sector. The 
gauge Bethe Ansatz equations read 

1 



^iLv, _Y[ ''^^Pk)-x (pj) 2x+{pk)x jpj) ^ 

j_fj^x-{pk) -x+{pj) ^ 

2x-{pk)x+{pj) 

Now, we consider L G 2N. 
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At g = the above equations reduce to those of the Heisenberg model. The Bethe 
momenta pk are non-trivial and can be determined numerically at each L as follows. At g = 
we have Uk = ^ cot The variables Uk can be determined by solving {e.g. iteratively) 

L 

27r Jfc = 2 arctan(M/j — uj) — 4L arctan(2 u^), (4-2) 
i=i 

where the Bethe quantum numbers {Jk} for the AF state are 

( L-1 L-3 L-3 L-l\ 
{Jk} = {-^'-^'•••'^'^1 • (4-3) 

Here, it is not useful to compute the weak coupling expansion of A. Indeed, due to the 
appearance of the non-trivial one- loop Bethe roots, all the coefficients of the expansion 
depend on L. 

On the other hand, we can integrate numerically the equations. The result is not 
surprising and could be expected on the basis of the Hubbard model solution at finite L 
discussed in [0]. At large A, all Bethe momenta flow to constant values pk Pk given by 

{p„...,P,} = [±-,±-,±-,...,±^.y (4.4) 
Hence, the asymptotic form of the anomalous dimension is 



A 1 \/A^p 

V 



. V 
sm — 
2 



(4.5) 



lim — = — + .... (4.7) 



^/A^-^ . p (2s + l)7r y/\ 

= > sm - = > sm = — 4.6 

Ltt ^ 2 Lit ^ 2L 27rXsin^ ^ ^ 

p>0 s=0 2i 

In particular, taking L ^ oo we find 

A_ _ ^/A 

L^'io 2L vr^ 

The factor 2L in the scaled anomalous dimension is the correct one, i.e. the length of the 
associated lattice model. Indeed, in the su(2) sector, the spin zero cyclic state associated 
with the AF state with L Bethe momenta has 2L spins, L with spin up and L with spin 
down. 

We remark that the above leading term is obtained both from the GBA equations and 
from the Hubbard model. It is quite interesting to see what happens at the level of string 
Bethe Ansatz equations. 

5. The string Bethe Ansatz equations for the AF state in the 5u(2) sector 

The string Bethe Ansatz equations in the su(2) sector are modified by the same universal 
dressing factor we introduced in the su(l|l) sector. Thus, they read 

1 



jLp, ^ TT a;+ {pk)-x jpj ) 2x+{pk)x jpj ) ^.^^^ ^ ^ ^^^^^ 

-'^ i — -. r —, r- 

2x~{j)k)x+{pj) 
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where the scattering phase ^{pk,Pj) has been defined in Eq. (|3.3D . 

We solve numerically these equations and the outcome is that aU Bethe momenta 
vanish hke pk a^A"^/^ precisely as in the su(l|l) case. Again we can work out an exact 
equation for the asymptotic coefficients {ofc}. Taking the limit of the SB A equation we 
find the following modified form of Eq. ( |3.11| ) 



exp 



an 



n 



a 



; (afc - aj) + 2-n{a^ ^ + a- 



a 



a 



h{a) = 7^ - 7^ 



,2 



Ai{h{ak)~h{aj)) 



a,- 



(5.2) 



a 

16tt' 



where in the l.h.s. the sum includes the case j = k. 

As in the 5u(l|l) sector, we can solve numerically this equation to cross check the 
numerical solution of the SBA equations. For instance, at L = 8 we have four positive 
vanishing momenta and the above equation predicts 

ai = 2.7192199579, (5.3) 
02 = 4.1578685742, 
as = 5.7295537708, 
04 = 12.5260989258. 

The actual comparison with the solution of the SBA equations is shown in Fig. (p!^). 
To find the asymptotic expression of the anomalous dimension we can follow the same 
strategy as we did in the su(l|l) sector. The extra factors in Eq. ( ^.2D also cancel when all 
the equations are multiplied together. Hence, we find again the fundamental relation 



E 



O-i 



(5.4) 



with a different sequence Nl- Evaluating the solution to Eq. ( p. 21) we find the table 

L : 4 6 8 10 12 14 16 18 20 ... 
iVi : 4 9 16 25 36 49 64 81 100 .. . 

Hence, the following simple formula holds 



(5.5) 



Nl 



(5.6) 
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leading to the prediction 

A string 

^-c,X^I\ (5.7) 
with ^ 1 L 1 

This is an exact result holding at any finite L. For instance, it can be checked that it is 
valid for the L = 8 solution up to the quoted accuracy. 

Again, the integer A'^^, appears in the asymptotic expression for A at large A in the 
form 

A~20vZA^/^ (5.9) 
identifying A'^^: with the level of the massive string state dual to the AF operator. 

6. Discussion 

Quantitative tests of AdS / CFT require a perturbative window allowing reliable calculations 
on both sides of the correspondence. Such a window does not exist for generic sectors of 
the spectrum, but is available for semiclassical string states with large quantum numbers. 
In these BMN-like limits, a perturbative check of AdS/CFT at weak effective 't Hooft 
coupling can be attempted, but is known to fail at three loops. This discrepancy can be 
seen as a limitation to our capability of capturing the strong quantum dynamics of string 
theory beyond the BMN limit, i.e. the small A regime for states with fixed finite classical 
dimension L. 

On one hand, the GBA equations are effective in computing all-loop perturbative gauge 
theory properties, like anomalous dimensions. However, the implicit order of limits (L — > oo 
after A — > 0) spoils the agreement with string calculations beyond two loops. On the other 
hand, the SEA equations are valid at large A, but already the leading quantum corrections 



are known to receive important and non-trivial corrections at small A [p3, 23, 24, 25], not 



yet under full control despite recent progresses |3^, 34 1. 

At large A, the picture is quite different. In principle, the GBA equations should not 
be trusted because the wrapping terms cannot be neglected. Instead, the SBA equations 
are expected to match string calculations, including leading quantum corrections. Hence, 
we can make predictions on the string side, but we cannot test the AdS/CFT correspon- 
dence. Actually, this general statement is very conservative. Specific cases must be treated 
with care and a matching between the solutions of the two sets of equations is possible. 
Indeed, the structural difference between GBA and SBA equations lies only in the dressing 
scattering phase Eq. (3^) and the relevance of this extra term should be considered case 



by case. 

Important examples of such exceptions can be found in various BA/[N-like limits. For 
instance, the exact all loop expression of plane wave string levels reads in the strict BMN 
limit 

^'^ I A 

Am- J = 5ZV^ + n2A', A' = -^, J^oo, (6.1) 

k=l 
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and scales like (A')^/^ at large A'. This result is valid both in the framework of GBA |11| 
and SBA [^] equations. In this case, the matching of the two predictions is due to the 
fact that the impurities are fixed in number and their diluteness prevents scattering effects 
in the thermodynamical limit. A related example is the Hofman-Maldacena limit which 
also displays asymptotic (A')-*^/^ scaling laws closely related to the BMN case [43, 44, 46 



45, 47|. Again, in the Hofman-Maldacena limit, the dressing factor in the SBA equations 
has been shown to decouple, making the prediction from the GBA exact, at least in the 



thermodynamical limit |48| 



These examples should be regarded as exceptions, precisely because the irrelevance of 
the dressing phase is not expected to be generic feature. A more involved example where 
an explicit strong coupling discrepancy appears is the folded string (FS) solution ||2^ in 
the 5u(2) sector. The energy of this solution is a function Afs(A, J). In the AdS/CFT 
correspondence, it must be matched with the anomalous dimension of an operator with 
L = 2J constituent scalar fields. The scaling operator is well known, at least in the L ^ oo 



limit and is the double contour solution of the GBA equations described in |37]. It has 



been studied in some details also in the Hubbard model formulation [14|. Setting, as usual 
A' = A/J^ we have 

Apg™^(A' J^, J) , 
hm ^— = /(A ), (6.2) 

j—foo ZJ 

where the function /(A') is explicitely known. The leading term is obtained by expanding 
at large A': 

/(A') ~ (A')^/l (6.3) 



One can ask whether it is possible to reproduce Eq. ( |6.3D with the GBA equations. The 
anomalous dimension A^g"^'' of the double contour solution is known at strong coupling 



and finite J in the Hubbard model GBA equations |14|. It reads 



A gauge / \ T\ ^ i 

With the usual remarks about the order of limits, we see that the GBA/SBA equations 
predict an asymptotic behaviour ~ (A')*^ with i' = 1/2 and 1/4 respectively. Here the 
number of impurities in the folded string solution gets large as J ^ cxd and their finite 
density makes the role of the dressing scattering phase non-trivial. 

Following these remarks, it seems interesting to look at other explicit non-trivial ex- 
amples where the GBA and SBA equations can be compared in the strong-coupling region. 
From a slightly different perspective, one is considering a special class of states and wonders 
about the role of the SBA scattering phase. Following this line of reasoning, this work an- 
alyzes the highest state in the su(l|l) and su(2) sectors of AA = 4 SYM. We have been able 
to solve all ambiguities appearing in the strong coupling expansion of the Bethe Ansatz 
equations. Our results have been cross-checked with a resummation technique which is 
able to connect smoothly the weak- and strong-coupling regions. Our main results have 
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already been summarized in Eqs. ( |l.5D and are repeated below for the reader's advantage. 

A™ 1 A string 

27rLsin— 2L 2 

(6.5) 



^string 

111 = AV2. 

L 

CL 0.1405(1), as L ^ oo. 



The main outcome of our analysis is the following. For any fixed L, the highest states 
in the su(2) or su(l|l) sectors have a large A anomalous dimension scaling like A ~ A'^ 
where u = 1/2 in the GBA equations and v = 1/4 in the SBA equations. At large A (and L) 
the SBA equations can be trusted without subtleties. Hence, the A^/^ scaling predicted by 
the GBA equations is not a true feature of the highest states. Their anomalous dimension 
immediately scales like X^/'^ and uniformly in L as soon as they are treated by the string 
Bethe Ansatz equations. 

We remark that our result is somewhat novel. Indeed, in the su(2) sector, it is common 
lore to believe in the ^/\ scaling of the AF operator, after its identification with the dual of 



the slow-string limit solution in [38|. From our analysis, we see that it is certainly possible 
to force the SBA equations to exhibit A^/^ scaling. However, this must be done by assuming 
a large A behavior of the Bethe momenta that is ruled out by the explicit solution of the 
equations, at least for the highest states. 

Notice also that it is possible to quantize the superstring equations of motion after 
truncation to the 5u(l|l) sector |4^, 50|. The spectrum contains long string solutions with 
non- vanishing winding w = ^^Pk with A^/^ scaling. On the contrary, short strings with 
vanishing winding exhibit the usual X^l'^ scaling. The observed symmetry p — > —p of Bethe 
momenta favors the w = ^ option. 

In conclusion, apart from the above mentioned special cases, it seems definitely dan- 
gerous to rely on the GBA equations to estimate the strong coupling limit of general states, 
as our analysis of the highest states has shown. Instead, the full solution of the SBA equa- 
tions, even at the discussed semi-analytical level, appears to be an effective predictive tool. 
For instance, our result ^ 

(6.6) 



string 



gives a simple formula for the level of the string state dual to the highest state in the two 
compact sectors. 



^"Nl is integer in both sectors since it has been derived with L odd (even) in su(l|l) (su(2)). 
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A. The Weniger resummation algorithm 



Given the power series 

A 



L—fOO L 

we can evaluate the partial sums 

Sn = >>fc- (A.2) 







n>0 












fe=0 




r approximants 




n\ (l+i)n-i 




j) (l + n)n-i 




n\ (l+j)n-l 


1 


j) (1 + n)n-\ 





(A.l) 



= ^ ' (A-3) 

where (a)m = r(fl + m)/r{a) is the Pochhammer symbol. If is beyond the convergence 
radius, the partial sums do not converge and oscillate wildly. For better stability we have 
performed all calculations in exact arithmetics. If the Weniger algorithms succeeds in 
resumming the series, then the Weniger approximants converge. 
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Figure 1: Wcnigcr convcrgcnts Snio^) for the weak coupling expansion of the GBA equations at 
the three values = 1, 10,30. 




Figure 2: Resummation of the weak coupling expansion of the GBA equations. We also show the 
result from the iterative numerical solution of the same equations at L = 15. 



L = 7 




Figure 4: Positive and vanishing Bethe momenta at L = 7. The dashed lines are the analytical 
predictions. 
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L = 9 




Figure 6: Positive and vanishing Bethe momenta at L = 9. The dashed lines are the analytical 
predictions. 
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L = 9 




Figure 7: A scaling of the vanishing momentum at L = 9. 
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L=ll 




Figure 8: Positive and vanishing Bethe momenta at i = 11. The dashed Unes are the analytical 
predictions. 
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Figure 9: Positive and (two) vanishing Bethe momenta at L = 13. The dashed hnes are the 
analytical predictions. 
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Figure 11: Determination of the prefactor Coo associated with the leading term in the strong 
coupling expansion of the anomalous dimension predicted by the GBA equations. As explained in 
the text, the three curves are subsequences of lattice length with {{L — l)/2) mod3 = 0, 1, 2. 
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Figure 12: Comparison of the gauge/string Bethe Ansatz solutions as predicted by the Weniger 
resummation algorithm. 
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Figure 13: Solution of the SBA equations at four L. We show a plot of X^^^p for the positive 
Bethe momenta. 
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Figure 14: Check of the finite L prediction of asymptotic coefficients ak from the solution of 
Eq. ( |3lll ). 
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Figure 15: Numerical clieck of the prediction Cao = l/\/2- 
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